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ABSTRACT 

Aims. We aim to perform the first long-term analysis of the system HS Hya. 

Methods. We performed an analysis of the long-term evolution of the light curves of the detached eclipsing system HS Hya. Collecting 
all available photometric data since its discovery, the light curves were analyzed with a special focus on the evolution of system's 
inclination. 

Results. We find that the system undergoes a rapid change of inclination. Since its discovery until today the system's inclination 
changed by more than 15°. The shape of the light curve changes, and now the eclipses are almost undetectable. The third distant 
component of the system is causing the precession of the close orbit, and the nodal period is about 631 yr. 

Conclusions. New precise observations are desperately needed, preferably this year, because the amplitude of variations is decreasing 
rapidly every year. We know only 10 such systems on the whole sky at present. 
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1. Introduction 

Eclipsing binaries are astronomical objects of high importance, 
especially owing to the possibility of deriving the basic physical 
properties of these stars with high precision. This is mostly be- 
cause one can easily calculate the individual masses, semimajor 
axis, etc. if one knows the inclination of the system and the ra- 
dial velocity curves. However, in some systems the plane of the 
orbit is moving slowly, and the radial velocity data have to be 
obtained at the same time as the data for the light curve solution. 
Otherwise, the method yields incorrect results. 

At present, we know only a few systems where the or- 
bital plane is moving and the eclipsing light curve had differ- 
ent shape s in dif f erent epochs. Six such systems were summa- 
rized by iMaverl (I2005I) . Moreover, three more systems show 
changes of minima depths, therefore they are also suspected 
to undergo a prec ession of the orbits, these are V685 Cen 
dMaver et al.L |2004|). AH Cep dDrechsel et all 1 1 9891). and V699 



Cvg (lAzimov & Zakirovlll99C 



2. The system HS Hya 

The eclip sing binary system HS Hya was discovered to be vari- 
able by Strohmeier et al.l (1 19651) . who also classified the sys- 
tem as an Alg ol-type, but the orbital period given is incorrect. 
Popped (119711) measured the radial velocities (hereafter RV) of 
the system, and analyzed the RV curves. The spectral type was 
derived to be F3-4, the correct orbital period is given to be about 
1.568024 days, and the mass ratio is about 0.96. However, the 
RV data were obtained over a period of five years, from 1966 
to 1970. Later, the complete light curve ( hereafter LC) was 
obtain ed in the Stromgren uvby system by iGvldenkerne et all 
(Il975h . These data were measured in 1972. The authors used 
the RV results from Pop per] (1 19711) and their inclination of about 
85.3°, which yielded a reliable picture of the system. 



Table 1. The inclination as obtained from various light curves. 



Year Inclination [deg] 



Reference 



1964 


88.9 


± 1.1 


Strohmeier et al. (1965) 


1972 


85.30 


+ 0.41 


Gyldenkerne et al. (1975) 


1991 


79.83 


± 0.21 


Perrvman etal. (1997) 


2002 


76.13 


±0.15 


ASAS 


2005 


75.19 


±0.28 


ASAS 


2008 


74.60 


±0.50 


ASAS 



However. iTorres et al.l (1 19971) published a new finding about 
HS Hya. The analys is was based on older data from photom- 
etry, RVs by iPoppen but also Torres and coworker's own new 
RV data, revealing that one more distant component is orbit- 
ing around the eclipsing pair. The period of this body is about 
190 days, and it is probably of the spectral type M0. Its light 
contribution is quite low (below 1%), but the RV residuals ob- 
tained by a cross-correlation clearly show periodic modulation. 

3. The change of inclination 

The star has also be en observed by the Hipparcos satellite 
dPerrvman et all IT997h . During its three-year mission both min- 
ima were observed, but the coverage is onl y poor. However, af- 
ter tra nsformation from H p to V magnitude (Harmanec & Bozic, 
2001), it is clear that the depths of both minima are much lower 
than in the LCs obtained 20 yr ago. We s olved the Hipparco s 
LC with the same para meters as given in | Tqrres et al.l dl 997b . 
The PHOEBE program dPrsa & Zwitteil 120051) was used, which 



is bas ed on the Wilson-Devinney code, Wilson & Devinnevl 
(fl97Tl) . All relevant parameters were fixed except for the incli- 
nation, see Table [TJ The value of the third light was fixed at a 
value of 0.4% only , in agreement with the finding published by 
iTorreset al.l(ll997l) . 
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Fig. 1. Available light curves of HS Hya in the V filter. The 
changing depth of both minima is clearly visible. The three bot- 
tom figures were plotted with the same range in the y-axis. 



The star was also in cluded into the photometric survey 
ASAS dPoimanskiL l2002h . We divided the whole data set into 
three parts and separately solved the light curves in 2002, 2005, 
and 2008. The procedure of LC fitting was the same as for the 
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Fig. 2. Inclination of the eclipsing binary with respect to the 
time. 



Hipparcos data, and the results are given in Table Q] As one can 
see from FigQ] the depths of the minima are still decreasing. 

Unfortunately, it is not easy to find other reliable photometry 
to do a similar analysis in different time epochs. One of the lim- 
iting problems is the role of the filter, because the entire above- 
mentioned photometry can easily be transformed into the stan- 
dard V magnitudes. Th e photometry fro m the automatic survey 
called "Pi of the sky" dBurd et all I2005I) is another possibility, 
but this photometry is unfiltered, and therefore it is problematic 
to solve its light curve. Moreover, it has fairly high scatter and 
covers a similar time span as the ASAS data. 

We also tried to us e the data from the discovery paper 
(IStrohmeier et all Il965l) . but these are only the photographic 
data and were not obtained in any standard photometric fil- 
ter. Another problem is that the original data are not avail- 
able, only the phase plot, but this was constructed with incor- 
rect ephemerides. As on e can see from the LCs published in 
iGvldenkerne et all d 1975b . the individual depths in different fil- 
ters are quite similar to each other. Moreover, fixing the other 
relevant parameters during the fitting process, we were able to 
construct a plot of minima d epth versus inclination. Using the 
eight dimmed data points from Strohm eier et al.l(ll965l) . we were 
able to roughly derive the inclination from these data obtained in 
1964, see TableQ] 

The data from Table Q] were used to construct the plot given 
in Fig|2] Fitting these data points with a linear curve, one sees 
that the change of inclination is about 0.3° during one year. 
Therefore, the amplitude of photometric variations is decreasing 
rapidly every year. 



4. The nodal period 

The precession effect of the close pair' s orbit due to the d istant 
third body was described elsewhere, e.g. ISoderhielmldT975l) . The 
nodal period can be computed from the equation 



nodal 



1 + 



Mi + M 2 



(1 



s!) 3 * (£ cos., | 



where subscripts 1 and 2 stand for the eclipsing binary compo- 
nents, while 3 stands for the third distant body. The term Go 
stands for the angular momentum of the wide orbit, and the 
C is the total angular momentum of the system. However, the 
problem of the unknown inclination of the wide orbit (which 
is included i n the last term in brackets) led us to use a differ- 
ent approach. lDrechsel et al.l(ll994l) analyzed the system IU Aur, 
where a similar problem arose, hence one can also fit the term 
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Table 3. The masses as derived from different data sets/methods. 
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Fig. 3. Fitting the cos(/) with the sinus term on available data, 
the confidence level of 95% is shown with the dash-dotted line. 

cos i with a sinusoidal fit, following the equation 

cos i = cos / • cos z"i - sin / • sin i\ ■ cos(2tt(/ - fo)/-P n odaiX 

where / is the inclination of the invariant plane against the ob- 
server's celestial plane, i is the inclination of the eclipsing bi- 
nary, while i\ is the inclination between the invariant plane and 
the orbital plane of the eclipsing binary. 

A similar analysis was performed (see Figj3), but regrettably 
only a small part of the nodal period is covered with data points 
nowadays. Fitting a sinusoidal curve, the nodal period is about 
631 yr, of which only about 1/14 is covered. The other two ad- 
justable quantities were only poorly constrained. New precise 
observations are urgently needed. 

5. Observational consequences 

HS Hya belongs to a group of unique systems, and therefore 
more observations are needed to confirm our hypothesis with 
higher confidence and to study its other physical properties. The 
amplitude of eclipses is lower than 0.015 mag in V filter at 
present, and is still decreasing. However, the ellipsoidal varia- 
tions outside the eclipse will remain even when the photomet- 
ric eclipses disappear. According to our model the photometric 
eclipses will stop in about 2022. Nevertheless, detecting these 
shallow eclipses is problematic, especially from the northern 
hemisphere, owing to the low declination of the star. 

Concerning the third body period, one can ask why the 190- 
day orbit was not discovered earlier via analyzing the minima 
times via period variations. There exists a huge database of min- 
ima observations (more than 100), but no variation was detected. 
This is because of the sh ort period of the third body and its low 
mass. This amplitude (see llrwinl 19591) resulted in about 1 minute 
only, which is comparable with the precision of individual times 
of minima observations. 

6. Discussion and conclusions 

Assembling all available systems with changing minima depths 
(see Table 13, one can see how unique these systems are. We 
know only 10 such systems today and the nodal period was de- 
rived in only four of them. Moreover, it seems that this effect 
was preferably observed in early type systems (B and A spectral 
types), and HS Hya is the first exception. 

Our hypothesis of changing inclination would also slightly 
shift the physical parameters of components as presented in 
iTorres etall (1 19971) . These authors assumed constant inclina- 
tion and used light curve and radial velocity data from different 





RV 1968 


RV 1992 


Torres etal.(1997) 


Mi/M e = 


1.319 


1.307 


1.255 


M 2 /M Q = 


1.291 


1.267 


1.219 



epochs. However, the time gap of more than 20 years between 
photometry and spectroscopy yields a difference in inclination of 
about 7.8°. This difference is able to shift the true stellar masses 
as computed from the term m-sin 3 (Q . The inclinat ion in 1968 
(when the RV data were obtained by Pop perl[l97ll) was about 
87.15°, while in 1992 (roug hly the middle o f the time interval 
of radial velocities used in ITorres et al. 1997) was about 79.3°. 
Using these values, the masses resulted in the values presented 
in Tabled 

These values agree much better with the predicted stellar 
masses of F4+F5 spectral types (i.e. with the temperatures) for 
norm al metallic i ty. He nce, the errors of masses as presented 
by ITorres et all (1 19971) of about 0.007 M Q are too optimistic, 
because they negl ect the systematic effect described above. 
ITorres et"ai] ([1997) presented a perfect fit of the derived values 
of log M, log T, and log g on the model isochrones. Likely an ad- 
justment of metallicity or age will be required to accommodate 
the new mass determinations. 

The plot presented in Fig. [3] needs to be spread in the next 
years. However, deriving the inclination of HS Hya when it stops 
having eclipses will be hard, because the ellipsoidal variations 
have only low amplitude. On the other hand, the interferometry 
of the close pair would solve this problem. Detecting the two 
eclipsing components via interferometry is difficult, but worth 
trying. The system itself is relatively bright and the two eclips- 
ing components have a similar luminosity (i.e. magnitude differ- 
ence close to zero). The computed angular separation of the two 
eclipsing components is about 0.4 mas. 
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Table 2. Known eclipsing binaries with changing minima depths. 



System 


Mag 


Sp. 


Eclipsing 


Long 


Nodal 


Reference 




V 


type 


period 


period 


period 






[mag] 




[day] 


[day] 


[yr] 




IU Aur 


8.39 


B3Vnne 


1.8115 


293.3 


330 


Ozdemir et al. (20031 


V685 Cen 


8.85 


AO 


1.1910 


7 


7 


Maver et al. (20041 


AH Cep 


6.88 


B0.5V 


1.7747 


7 


7 


Drechsel et al. (19891 


V699 Cyg 


11.6 


B2 


1.5515 


? 


7 


Azimov & Zakirov (19911 


SVGem 


10.57 


B3 


4.0061 


? 


7 


Guilbault et al. (20011 


HS Hya 


8.07 


F5V 


1.5680 


190 


631 


This study 


SS Lac 


10.12 


B9V 


14.4162 


679 


600 


Torres (20011 


AYMus 


10.35 


B9 


3.2055 


? 


? 


Soderhielm (19751 


RWPer 


9.68 


A5Ve 


13.1989 


7 


7 


Olson etal. (19921 


V907 Sco 


8.61 


B9V 


3.7763 


99.3 


68 


Lacv et al. (1999) 



4 



